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We consider linear arrays of cells of volume 14 populated by monodisperse rods of size crl4, 
cr = 1,2,..., subject to hardcore exclusion interaction. Each rod experiences a position-dependent 
external potential. In one application we also examine effects of contact forces between rods. We 
employ two distinct methods of exact analysis with complementary strengths and different limits of 
spatial resolution to calculate profiles of pressure and density on mesoscopic and microscopic length 
scales at thermal equilibrium. One method uses density functionals and the other statistically 
interacting vacancy particles. The applications worked out include gravity, power-law traps, and 
hard walls. We identify oscillations in the profiles on a microscopic length scale and show how they 
are systematically averaged out on a well-defined mesoscopic length scale to establish full consistency 
between the two approaches. The continuum limit, realized as 14 —>■ 0, cr —>■ (Xj at nonzero and finite 
cr 14, connects our highest-resolution results with known exact results for monodisperse rods in a 
continuum. We also compare the pressure profiles obtained from density functionals with the average 
microscopic pressure profiles derived from the pair distribution function. 


I. INTRODUCTION 

In classical statistical mechanics, particles with shapes 
are ubiquitous. Their prominence in granular matter 
mm, soft condensed matter Ellll, and, more specifically, 
biological matter mm, is well established. Their shapes 
vary from the highly complex such as folded proteins to 
the most elementary such as hard spheres. Their equilib¬ 
rium and nonequilibrium properties are investigated by a 
broad array of experimental and computational probes. 

Analytic approaches in this area of research see their 
predictive power restricted to fairly simple scenarios re¬ 
garding shapes, environment, and dimensionality. For 
rigorous calculations the domains of applicability are fur¬ 
ther narrowed. This limitation finds ample compensation 
in their usefulness as benchmarks and anchor points for 
approximations and simulations. The work reported in 
the following is motivated by this chain of reasoning. It 
deals with hard rods in one dimension at thermal equilib¬ 
rium in external potentials. We consider monodisperse 
rods of size aVc, u = 1, 2,... on a lattice (linear array of 
cells with volume 14). Hard rods of size 14 populating a 
continuum emerge from the limit 14 —>■ 0, tr —>■ oo with 
(tUc = Ur. 

There are several approaches that facilitate an ex¬ 
act derivation of thermodynamic and structural proper¬ 
ties for a homogeneous one-dimensional hard rod system 
(Tonks gas) with first-neighbor Takahashi-type interac¬ 
tions [7]. One rather elegant method uses convolution 
relations between Boltzmann factors to determine parti¬ 
tion functions [5]. It produces the free enthalpy and the 
equation of state (EOS) for first-neighbor interactions of 
arbitrary range. Many-body-density distribution func¬ 
tions can be calculated by an extension of this approach 
[5]. An alternative method of similar scope has recently 
been developed. It uses statistically interacting vacancy 
particles (SIVP) as quasiparticles [10]. This method also 
yields the size distribution of vacancies between rods. 


On the basis of these approaches for homogeneous sys¬ 
tems, it is possible to treat inhomogeneous systems by 
assuming that the EOS is valid on a coarse-grained local 
scale. From the requirement of mechanical equilibrium, 
spatial variations of pressure and density can then be 
calculated from the spatial variations of the external po¬ 
tential. This provides a simple and common thermody¬ 
namic route for calculating density and pressure profiles 
m, which generally is a difficult task for an interacting 
many-particle system. We will refer to this thermody¬ 
namic route as the EOS method in the following. Be¬ 
cause the EOS method relies on the assumption of the 
existence of a local EOS, it is interesting to gain insight 
into how far this assumption is justified and whether the 
method gives useful information even if the underlying 
assumption does not hold. 

To tackle this question analytically, exact results for 
inhomogeneous systems are required. For hard rods with 
first-neighbor interactions in one dimension, exact treat¬ 
ments are possible via recursion relations for partition 
functions or density functional theory (DFT) [T3|- 

HSj. These methods allow for the exact derivation of 
density profiles as well as many-body distribution func¬ 
tions. Given exact density profiles, pressure profiles can 
be obtained by resorting to the requirement of mechan¬ 
ical equilibrium as in the EOS method, but without as¬ 
suming a local EOS. 

The exact calculation of local pressures is a more subtle 
task. Generally, the local pressure can be defined via the 
trace of the local stress tensor, which governs the time 
evolution of the momentum density in a coarse-grained 
continuum description mini. In case of pair interaction 
forces, thermodynamic averaging over the corresponding 
local stress tensor allows for the determination of pres¬ 
sure profiles from the density profiles and pair distribu¬ 
tion function. 

Our calculations here will employ two approaches. The 
DFT for lattice fluids jlH [TU is used to determine 
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exact density profiles in external fields and the SIVP ap¬ 
proach [ini ESHlTj is used as a realization of the EOS 
method. The DFT and SIVP approaches have domains 
that partially overlap and strengths that complement 
each other. In the DFT, the operational degrees of free¬ 
dom are the rods themselves. In the SIVP approach, the 
operational degrees of freedom are the vacancies between 
the rods. 

We begin by describing the general methodology and 
background (Sec. 0 and then proceed with appl ications 
to rods in a uniform gravitational held (Sec. Ill I and in 
a power-law trap (Sec. |IV[). The subtleties regarding av¬ 
erage microscopic pressure are addressed in the context 
of the Hrst application. Steric wall effects in lattice sys¬ 
tems and their relations to known continuum results are 
discussed in Sec. m In Sec. IVII we summarize the main 
conclusions and outline projected extensions to polydis- 
perse rods. Appendices and summarize outlying 
background materials for use in the main text. Appendix 
[^presents a highly practical method of calculating exact 
density prohles within the DFT framework for arbitrary 
external potentials. 


II. METHODOLOGY 
A. Model system 


The model system analyzed in this work is specihed by 
the Hamiltonian, 


'H(n) = ^ VijUiUj 

i<j 


lAiJii , 


(1) 


where Hi is an external potential and Vij an interac¬ 
tion potential restricted to hrst-neighbor rods. The hrst- 
neighbor restriction is naturally ascertained by interac¬ 
tions Vij with a range limited to j = f -I- where 
cr < ^ < 2(7, but an extended range between hrst- 
neighbor rods is permitted. 

In this work we mainly examine the effects of hardcore 
repulsion in combination with external potentials. Con¬ 
tact forces are included in one application (Secs. HID 
and HIE). The effects of long-range forces are being an¬ 


alyzed in a separate study [28] . 


B. Exact density functionals 

The analysis carried out in |2Ij is based on former work 
[n mulls and expresses the grand potential as a den¬ 
sity functional, i.e. a functional of the mean occupation 
numbers of rods hi = (rii) |3S 

L 

D[hi,..., hi] = F[hi,..., hi] -I- '^{Ui - , (2) 

i=l 


Consider rods on a linear chain i = ,L of lattice 
sites, represented in Fig. as an array of cells. Each 
cell has volume 14. Rods of size a occupy that many 
adjacent cells. Hard walls at both ends of the chain dehne 
the boundary conditions. Assuming that cells (and rods) 
have unit cross section we can conveniently use 14 as a 
microscopic scale for both volume and length. 

Microstates of this system are encoded in a sequence 
of occupation numbers, n = {ni,..., nz,}, rii = 0,1. To 
a rod that occupies sites i,... ,i -\- a — we assign the 
occupation number rii = I. Hardcore exclusion imposes 


the conditions 


= 0, j = 1,..., 

(7-1, 

and the 

bard walls imply 
respectively. 

rij = 

0 for i < 1 and 

i > L 

— fJ + 1, 
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FIG. 1: Four rods of size a — 3 a,t positions i = 2,7, 14,17 
on a lattice of size L — 20. The relevant external potentials 
Ui and interaction potentials Vij are stated. One is a contact 
potential and the other represents an interaction of maximum 
range (a - 1)14- 


where /i is the chemical potential. The intrinsic free- 
energy functional in Eq. ([^ can be written in the form 

i=l i=l 

where T is the temperature, and /i[h] = ei[n] — rsi[n], 
ei[n], and Si[h] are local functionals of the intrinsic free 
energy, internal interaction energy, and entropy. The lat¬ 
ter are given by 

i+i 

ei[n] = (4) 

j=i+<y 

and 

Si[n] = - fcB|4>(ai[n]) + $(&i.i[n]) - [n])-h 

i—(T 

+ [®(C'ij[n]) + 4'(d*j[n])-4>(d*-ij[n])]| (5) 

3=i-i 

with pair correlators 

Qj[n] = (niTij), (6) 

$(x) = x\a.x, and |37| 

i—(7 

a4n]=hi- ^ Cij[vL], 


(7a) 
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i i k—(7 

^*j[n] = l- nk+ (7b) 

k=j-i k=j-(+a- l=j-^ 


d^A^]='nj- ^k,j[n]- (7c) 

k=j+a 

The pair correlators Q have their dependence on n en¬ 
coded in the implicit relations, 




r< — 

— 7 I I J u 

’ k—i-\-l ’ 


n 


dk,jbk-l^k g-/3Vi , 


P = 


k^T' 


( 8 ) 


The equilibrium density profile of rods follows from the 
extremum condition, 

^f^[n] = 0, i = (9) 

With the solution n of (|^ , the functions fi [n], [n] and 

Si [h] become the intrinsic free energy, internal interaction 
energy, and entropy per site. The equilibrium density 
profile n and the profiles of the thermodynamic poten¬ 
tials depend on temperature, crowding, interaction, and 
environment via /3, n, Vij, and Ui, respectively. The 
cell occupancy (mass density) is obtained from the rod 
occupancy (number density) via 


Pi = 

j—i—(7+1 

In homogeneous systems the pressure p follows rig- 
orously from the free-energy density / = F/LVc via 
p = ndf/dn — f. A natural extension of this relation 
to systems with external potential has the form 


i 

Pi = -fi[h]+ 

j=i—(7+1 


dfijn] 
' dfii 


( 11 ) 


and produces pressure profiles on a microscopic length 
scale. However, there is no guarantee that the pressure 
thus derived coincides with the average microscopic pres¬ 
sure as commonly defined via the pair-distribution func¬ 
tion. More on this question follows in Secs. |IIE 

Exact profiles for hi and pi for arbitrary external po¬ 
tentials and interactions limited to hardcore repulsion on 
the lattice are calculated by the method introduced in 
Appendix [B| 


and III C 


C. Coarse graining 


For this comparative study of methods we need a con¬ 
tinuum description for rod positions on a micros copic 
length scale used in the DFT approach (Sec. IIBI that 
carries over naturally to the mesoscopic length scale used 


in any of t he EO S methods, specifically the SIVP ap¬ 
proach (Sec. IID). This continuum description of the lat¬ 
tice system is unrelated to the continuum limit. We re¬ 
place each lattice site i with the interval [iI4, (* + l)b(;[ 
across one cell (of unit cross section), and we define the 
interaction potential v(x,x') and external potential u{x) 
by setting v{xi,Xj) = Vij and u(xi) = +. The first- 
neighbor restriction for the interaction potential becomes 
v{x,x') = 0 for \x — x’\ > aVc- The local coverage at 
equilibrium of this interval allows us to define a number 
density by the piecewise continuous function 


1 ^ 

p{x) = — = - W hj, X e[iVc,{i + l)Vc[. (12) 

a G 

2—i—(7+l 


For a homogeneous situation we have p{x) = n. 

From any density profile p+) we can calculate the as¬ 
sociated pressure profiles p{x) by invoking the balance 
between internal and external forces at thermodynamic 
equilibrium 

ydp^ = p{x)fuix ), (13) 

da; 

where + = —du/dx is the external force field [5B]. Inte¬ 
gration of this differential equation yields 

p{x) - p{xo) = ^ [ dx'p{x')fn{x'). (14) 

JxQ 

In some applications the reference pressure p{xo) is 
known, e.g. via the weight of the rods in a uniform grav¬ 
itational field. In other cases it can be determined from 
the average number N of rods, which we know from sum¬ 
ming fij over all sites, by using the normalization relation 



dx p{x) 


, p{x) dp 

dx-r-r^, p = 

/u(x) dx 


(15) 


If xq = LVc then we can use p{xo) = —9E[n]/9a;o. If 
the rods are only subject to hardcore repulsion we can 
use the fact that at the system boundaries we have kine¬ 
matic pressure k+Tp (Sec. |IIE ). In Sec. IIIB we compare 
profiles inferred from 0 with profiles calculated from 
Eq. pi). 


D. EOS method and SIVP approach 


If one assumes that the EOS of a homogeneous system 
remains valid in a corresponding inhomogeneous system 
on a coarse-grained local scale in the presence of exter¬ 
nal potentials, the balance equation (13) is sufficient to 


determine density and pressure profiles. Depending on 
the circumstances we use (13) to calculate the functions 
p{x) and p(p(x)) or the functions p{x) and p(p(a;)) |39j . 

In the former case we have an EOS in the form p{p) 
and solve by separating variables p and x: 


Vc [ fAx') = u{xo) - u{x) 

Jpo pip ) Jxo 


(16) 
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where po = p(xo) is determined by one of the conditions 
stated in Sec. |II C| In the latter case we proceed analo¬ 
gously via separation of variables p and x. 

The EOS method is particularly useful if long-range in¬ 
teractions are present. For such cases an exact DFT cal¬ 
culation of density profiles tends to be be impracticable 
and calculations based on recursion relation for partition 
functions are cumbersome. For long-range first- 

neighbor interactions, the SIVP method m provides a 
user-friendly way to derive the EOS. 

The microstates are encoded in a sequence of — 1 
vacancies of size m (m = 0,1,2,...) between consecutive 
rods, N in number. Summing over all microstates means 
summing over all size combinations of iV — 1 vacancies. 
This sum is free of constraints. The interaction energy 
of first-neighbor rods at distance m is equivalent to part 
of the excitation energy e-mip) of the vacancies, 

^m{p) = pm + (j)m, (17) 


where 4>m = Vi.i+m+o-- The vacancies themselves are 
free of interaction energies and form a set of polydisperse 
quasiparticles with generalized exclusion statistics. Their 
statistical mechanics has been worked out in m building 
on a host of foundational work [52H1Z]- This treatment 
produces exact results for any thermodynamic quantity 
of interest for spatially homogeneous situations. 

The free enthalpy G{p) per site is given by 


/3G(p) 


1" l+F 




(18) 


from which other thermodynamic quantities are inferred 
via the auxiliary quantities |10j 


Baj{p) = m°‘[/3em{p)re (19) 

m—O 


This includes the mean size m of vacant cells 

mW = ' 

Boo[P) 

from which the EOS is obtained in the form 

1 


Kp) = 


m(p) 


( 20 ) 


( 21 ) 


E. Local pressure from pair density function 


In continuum mechanics, the ‘microscopic’ pressure is 
defined as one third of the (negative) trace of the mi¬ 
croscopic stress tensor and the divergence of this tensor 
determines the time evolution of the momentum den¬ 
sity caused by the internal interaction forces. In one- 
dimension, the local stress tensor reduces to the micro¬ 
scopic pressure Pniic(x,t) and the equation of motion for 
the momentum density n(a;,f) becomes 


_ ^Pmic „ext 

dt ~ dx ^ 


( 22 ) 


How statistical mechanical expressions for the local pres¬ 
sure Pmicix) = (pmic(a:,t))eq in thermodynamic equilib¬ 
rium (or the equilibrium-averaged microscopic stress ten¬ 
sor) are obtained when starting with these Euler equa¬ 
tions of continuum mechanics was first studied by Irv¬ 
ing and Kirkwood in 1950 [16]. In Appen dix [A| we have 
adapted the elegant derivation by Lutsko |17j to one di¬ 
mension, which for pair interactions yields 

Pmic{x) = ksTp{x) + dxi dX2P^'^Hxi,X2)f{xi,X2) ■ 

Jo Jx 

(23) 


Here p^^'>{xi,X 2 ) is the pair distribution function and 
f{xi,X 2 ) the force of a particle at position xi on a parti¬ 
cle at position X 2 - The first term represents the kinematic 
pressure and the second the interaction pressure. We set 


He = 1 throughout Sec. HE 


For hard rod systems, the interaction potential v = 
v(\x 2 — Xi\) is infinite for \x 2 — Xi\ < a. This singularity 
can be accounted for in Eq. (23) by considering a mod¬ 
ified continuous potential v^yP}, which agrees with v{r) 
for r > a while for r < ct it is given by Ve{r) = vq for 
r < a — e and v^{r) = v{a) — [('Cq — 'v{a))/f\{r — a) for 
a — e < r < a. After inserting the corresponding force 
in Eq. (23) the pressure is obtained by taking the limit 
—>■ oo, e —)■ 0 . For non-interacting hard rods (n(r) = 0 
for r > (t) in particular, this procedure yields 

Pmicix) = kBTp{x) + kBT f dx'p^'^'> (x', x'+ cr). ( 24 ) 

J X — G 


Because (a;, y) = 0 for a; < 0 or ?/ > L — cr, the range 
of integration in the second term extends from zero to x 
for X < a, and from a; — cr to L — 2cr for (L — 2a) < x < 
(L—a). This means that the interaction pressure at a; = 0 
and X = L — a (the effective system boundary) vanishes. 
Accordingly, the local pressure at these boundary points 
is just given by the kinematic pressure, as earlier pointed 
out by Ibsen et al. [29|, who derived Eqs. (23) and (24) 
based on the approach in m- For interacting hard rods 
(H(r) 7 ^ 0 for r > a), the interaction pressure generally 
does not vanish at the system boundaries. 

It is interesting to note that equating the expressions 
for the local pressure in Eqs. (14) and (23) [or Eq. (24)] 
yields an integral equation connecting the pair distribu¬ 
tion with the density. This could in principle be used to 
determine p^‘^\x,y). Alternatively, the pair distribution 
can be obtained by solving the inhomogeneous Ornstein- 
Zernike relation with the direct correlation function given 
by the second order derivatives of the density functional. 
For hard rods with first-neighbor interactions, many- 
particle-density distribution functions are most conve¬ 
niently obtained by employing recursion relations for par¬ 
tition functions, both for continuum |9| and for lattice 
systems [12]. We note that for the latter, the integrals 
for the interaction pressure in Eqs. (23) an d (|M| can be 
replaced by corresponding sums (Appendix |C[). 
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III. GRAVITATIONAL FIELD 


Consider a semi-infinite vertical column of cells num¬ 
bered 1 = 1,2,... from the bottom up. A uniform grav¬ 
itational field g acts on rods of mass m^. We use it 
here to represent any linear potential. We begin with 
non-interacting rods of size cr on a lattice and then take 
the continuum limit. Results from SIVP operating on a 
mesoscopic length scale are compared with those from 
DFT operating on a microscopic length scale. Pres¬ 
sure profiles on a microscopic length scale obtained from 
DFT via (11) are then compared with average micro¬ 


scopic pressure profiles inferred via (24). Finally, we dis¬ 


cuss some effects of repulsive contact interaction as made 
manifest in one or the other method. For the sake of 
brevity we limit the discussion to one simple case study 
of each approach. They can both be adapted to different 
applications. 


A. SIVP approach 


The gravitational potential to be used in (13) is 
u{z) = m^gz. 


(25) 


Convenient scaled variables for position, pressure, and 
temperature in this application are 

z=—, P=—, T= -— , (26) 

Ps PsCrVc 

where Zs = NaVc is the length of all N rods stacked up 
in a solid column. The pressure at 2 = 0 is Ps = Nm^g, 
independent of T. The thermal energy is measured 
in units of the work PsCrVc required to lift this weight a 
distance equal to the size o f one rod. 

The SIVP analysis (Sec. [IIP I starts from the expres¬ 
sion for the density of vacant cells, 


m = 


exp 


^ 1-1 
aT 


(27) 


derived in m from ( (20| . The scaled mass density (vol¬ 
ume fraction) inferred from (21) then reads 


=ap = 


(*) 


exp - 1 


exp 


(*) 


-^ + h 


(28) 


Performing the integral (16) with (25) and (281 yields the 


following equation for the pressure profile 


crT 


o1/<tT _ I 


(29) 


The density profile follows from (28) by substitution. 

The solid curves in Fig. show the profiles p{z) and 
p(mes)|'£^ at various values of T for rods of size a = 1. 



z 


z 


FIG. 2: (Color online) Profiles of (a) pressure and (b) density 
of rods in a uniform gravitational field at different tempera¬ 
tures. The solid lines represent rods of size cr = 1 on a lattice 
and the dashed lines rods (of any size) in a continuum. 


The variation of pressure with height crosses over from 
hydrostatic to atmospheric with increasing temperature. 
The density profile is rectangular in the low-temperature 
limit and varies like pi™®®) -wp/TatT^l. These pro¬ 
files do not vary much with the rod size a when expressed 
by the scaled quantities ( [2^ . In the limit cr —oo, 14 —f 0 
with crI4 = 14 fixed, we have a system of rods of size Vi¬ 
and mass iTir in a continuum. The continuum versions of 


(28) and (29) read 


(mes) _ 


1 


T/p+l 


p — 1 -\- Tlnp = —z, 


(30) 


respectively. The pressure and mass density profiles in 
the continuum are shown as dashed curves in Fig. 


B. DFT approach 

On the length scale of single rods additional features, 
not resolved by any EOS method including SIVP, emerge 
in the profiles for cr > 2 when analyzed via DFT. We 
write 


Ui = rurgzi 


Zi= 


(31) 


where Zi is the position of the center of the lowest cell 
occupied by a rod. The free energy functional ([^ with no 
interaction except hardcore repulsion acquires the form 
(Bl) and the density profile of rods for any given Ui are 


the solutions of the coupled equations (B21. Here, for the 
linear potential we set M —> 00 in all expressions 
imported from Appendix [B| 


We first examine the case cr = 1. The solution (B4) is 


constructed from an exponential function as follows: 


Ui = 


CA* 


1 -I- CAi 


A, = 




(32) 


where Zi = Zi/{N)aVc and the average number (N) 
of rods is controlled by the fugacity 4 = where 

fi = p/psCrVc is the scaled chemical potential. This DFT 
result exactly reproduces the functional dependence of 































6 


p(mes) - obtained via SIVP and given by (28) with 
(29) for cr = 1 if we use (10) and set 


c = - 1. 


(33) 


The results from both methods are fully consistent. The 
SIVP solution remains exact even for small numbers of 
rods. There exist no microscopic features in the density 
profile that SIVP does not resolve. 

Now we turn to the case <7 = 2, where the microscopic 
length scale does indeed reveal additional structures in 
the v ario us profiles. These structures are encoded in 
Eqs. (B5), which for slowly varying profiles we expand 


into the form, 


CA.= 


njjl - Hi) 
(l-2h,)2 


1 + 


rii-i — 2ni + 

1 — 2fi,- 


, (34) 


and note that the first correction is of second order. The 
leading term alone leads to the density profile, 


pf= 2fi, = 1 - 


1 


which coincides with 


V1 + 4CA.’ 

and for cr = 2 if we set 


(35) 


_ l/2f( l/2f _ 


C = e 


!)■ 


(36) 


Thus full consistency between the two approaches is es¬ 
tablished on the mesoscopic length scale. 

Finding the microscopic structures in the density pro¬ 
files of rods and mass requires that we solve Eqs. (B5) 
and then use ( [To| ) instead of Eq. (35). The solution as 
derived in Appendix [B] reads 


OO 

E' 

1=0 


_iy TT 


z = I,2,..., 


with the hi determined recursively from 


— 0 , hi — 


CA* 


I + hi-i ’ 


i = I,2,.... 


(37) 


(38) 


The circles in the panels on the left in Fig.|^are derived 
from (371. The probability hi that a rod occupies cells 


i and i + 1 varies with index i in a manner that reflects 
the combined effects of the hardcore exclusion interaction 
between rods and the presence of a hard floor i = 0. The 
spatial oscillations are mild at high T and discernible only 
very close to the floor. As T is lowered, the amplitude 
becomes stronger and the range wider. In the limit T —>■ 0 
the hi strictly alternate between one and zero, reflecting 
a compact stack of rods. 


When transcribed via (10) to the mass density, the wall 


effect is not nearly as strong. The evidence is represented 
by the circles in the panels on the right, where we set 


z = +Zi-i) for the discrete data and z from (26) for 


the curves. At high T the effect is still strongest m the 
immediate vicinity of the floor but that is no longer the 



FIG. 3: (Color online) Density and pressure profiles for rods 
of size CT = 2 in a uniform gravitational field at differ ent tem¬ 
peratures. The solid curves represent solutions of (281 and 


([29| as predicted by SIVP for a system with N ^ 1. The hi- 
data on the left (circles) originate from ( |37[ ) with ( from | |36[ ) 
and Ai from ( |32[ ) with (N) = 5. These data are transcribed 
into the circl es o n the right via (101 and into the circles of 
panel (g) via (391. 


case at the lowest T used in Fig. At T = 0 the effect 
disappears altogether. 

The local pressure predicted by DFT as inferred from 
0 with the Hi from ( |B2[ ) substituted into the free- 
energy functional (BI) becomes 


PiVc = ksTln 



(39) 


with Pi from ( |I0[ ). 

In panel (g) of Fig. we show the pressure profiles 
thus obtained from the density data of panels (a)-(f). 
The profiles are monotonically decreasing from p(0) = I. 
Remarkably, the wild oscillations of the hi are almost 
completely smoothed by (39). However, there does exist 
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a systematic albeit small deviation between the DFT and 
SIVP pressure profiles. The DFT pressure profile is closer 
to but not identical with the true microscopic pressure 
as will be further discussed in Sec. IIII Cl 

The agreement of the DFT and SIVP approaches on 
the mesoscopic length scale defined earlier is underlined 
by expressions (|32[ ) and (33) for cr = 1 and by expres¬ 
sions (35) and (^36| for cr = 2. It can further be shown 


that for any a the density determined by (28) 

and ( |29[ ) is the solution of a polynomial equation of or¬ 
der a. Likewise, the probabilities fii at Zi inferred from 


(B2) with all hj within the square brackets set equal to 
each other is also the solution of a polynomial equation 
of order a. The control variables are T, iV in the first 
polynomial equation (canonical ensemble) whereas they 
are T and C in the second polynomial equation (grand- 
canonical ensemble). The two polynomial equations are 
equivalent if we set 


(mes) _ 

Pl = CTUi 


= p('r-l)/<rf / l/^rf _ 


C = e 


I). (40) 


The agreement between DFT and SIVP on the meso¬ 
scopic length scale also extends to the pressure profiles. 
From (39) with = afii we infer 


P 


’mes) 


—— In 1 + . ^ 

crT V cr(I-p(“®®)) 


(41) 


which is equivalent to (28). 


The microscopic features in the density profile remain 
conspicuous for cr > 2. The oscillations that are superim¬ 
posed on profile predicted by SIVP are characterized by 
a ‘wavelength’ proportional to cr. Such profiles are read¬ 
ily produced from (B15) and (B16). The exact hard-wall 


effects for rods of arbitrary size cr will be investigated 
in Sec. |V] Soft walls as realized in power-law traps also 
produce structures on a microscopic length scale. Some 
examples will be investigated in Sec. |IV[ 


C. Average local microscopic pressure 

The average local microscopic pressure Pmicix) as in¬ 
ferred via (24) from the pair distribution function by the 
method presented in Appendix [C| has two parts: a kine¬ 
matic pressure and an interaction pressure (in a formal 
sense). In Fig. we show the profiles of both parts 
and their sum for rods of size ct = 10 at low temper¬ 
ature (T = 0.1). Both parts show strong oscillations 
that are somewhat out of phase. These oscillations are 
strongly attenuated with distance from the hard floor. 
They quickly become imperceptible with rising T (see 
Fig. §. The non-monotonic features of Pmicix) remain 
totally unresolved in the EOS pressure profiles discussed 
earlier but are partially resolved by the DFT pressure 
profiles as illustrated in Fig.[^ At high T the two profiles 
are virtually identical and equal to the SIVP profile. At 
low T, the additional pmicix) oscillations are reproduced 
by DFT with remarkable accuracy albeit not exactly m- 



FIG. 4: Profiles of kinematic, interaction and total pressure 
Pmic for rods of size cr = 10 at low T as inferred from wM I and 
the pair distribution function as derived in Appendix^ 



O 

II 


o 

II 

0.8 

\ (7= 10 

0.8 

A* (7=10 


c 

II 


s 

II 

o 

0.6 


0.6 


0.4 


0.4 


0.2 


0.2 






“c 

1 2 3 4 ( 

0.5 1.5 


FIG. 5: (G olor online) Profiles of th e lo cal pressure calcu¬ 
lated from (241 in conjunction with (G7l (circles connected 
by dashed lines) and calculated from (solid lines) for 
high T (left) and low T (right). 


D. Contact interaction via SIVP 

Consider a contact interaction potential v that is at¬ 
tractive for u > 0 and repulsive for v < 0. As inferred 
from m, the contact interaction changes the local den¬ 
sity of vacant cells from (27) to 


1 


m = 


(1 - tP)[l -b (1 - tP)t-P-^] 
where we have introduced the variable 

t = (0 < t < 1) 


(42) 


(43) 


and the scaled in tera ction v = vjapsVc in addition to the 
scaled variables (26). The mass density is inferred 


from (21). The integral (16) can still be evaluated exactly 
and yields the pressure profile. 


t 


1 - t(i - r) 1 - tp 
1 - tp(l - r) l-t 


= t° 


(44) 


in generalization of (29). An attractive contact interac¬ 


tion affects the profiles in a way similar to what a drop in 
temperature does (see Fig. [^. No significant additional 
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FIG. 6 : (Color online) Profiles of (a) pressure and (b) density 
of rods (of size cr = 1 ) in a uniform gravitational field at 
temperatures T = 0.05 and for (repulsive) contact interaction 
D = 0, —0.5, —1, —5. Panel (c) shows the p versus z data near 
z = 0. The dashed line has slope —2/3. 


features make their appearance. Repulsion is more inter¬ 
esting in that respect as documented in Fig.|^for a = 1. 

The configuration of rods becomes stratified at low 
temperature. The density profile now exhibits an ad¬ 
ditional layer of intermediate density. The width of that 
layer increases with the strength of the repulsion. The 
(largely hydrostatic) pressure profile acquires different 
slopes inside different layers. 

These profiles depend systematically on a without pro¬ 
ducing any additional features. From m we know that 
the effects of attractive or repulsive contact interactions 
of finite strength fade away completely in the continuum 
limit. 


One subtle feature of note concerns the initial slope of 
the hydrostatic pressure as presented in panel (a) of Fig. 
and then again in panel (c) on a much expanded scale. 
What in panel (a) looks like a clean change between slope 
— 1 for interaction strengths D > — 1 and slope — 1/2 for 
stronger contact repulsion, h < — 1 is, at very small z, a 
two-step change with slope —2/3 at D = — 1 in the middle. 
This asymptotic slope is evident in the expansion of (441: 


p = 1 — z 



1 


0{z^). (45) 


E. Contact interaction via DFT 


Here we wish to demonstrate how a repulsive con¬ 
tact force of infinite strength effectively increases the 
size of rods. The inclusion of a scaled contact potential 


Vc = —v/T generalizes (Bl) to 

^^ = e| (i- E E 

i=l [ \ j=i-o- J \ j=i-<T 

Pi ll^ (^2 4“ '^i—a ln {j^i — a 

Inhi-o- - 1 - E ^ M “ E 


J=l-rT 


j=t-a 



as shown in | 21 j , where the contact interaction Vc is con¬ 
tained in the correlators. 


^ \/K - 4e ’'Ke "<= - 

2(e-"<^ - 1) ’ 


Ai — 1 -|- e {ui—cr fli) — 


k—i—cr 




(48) 


that appear in the first three terms of (46 1 . The last two 
terms reflect the hardcore repulsion. 

The metamorphosis of hard rods of one size into hard 
rods of a bigger size is most transparent if we consider 
the case cr = 1 and compare the limits Uc = 0 and Vc = 
4 - 00 . The extremum condition © applied to leads to 
the following set of relations that determine the density 
profile {hi} for external potential Ui and fugacity 


Ce 


-l3Ui _ 


(49) 


[1 - ni\[ni - - Ci^i+i] 


41 - rii-i -ni + - Ui- Ui+i + Cj^i+i] ’ 


^ _ A^- -477(77-4 l)n,_ini 

1,2 — 

277 

where Ai = l + rj(ni-i +ni) and 77 = — 1. For Uc —>■ 0 

we have Ci^i+i = hihi+i and (49) reduces to 


Ce 


-/3Wi _ 


(51) 


whereas for Uc —>■ 00 we have Ci^i+i = 0 and (49) becomes 


Ce 


-/3Ui _ 


nt(l - rii) 


(1 - n^-i - ni){l -rii- n^+i) ’ 


(52) 


representing rods of size cr = 2 with only hardcore repul¬ 
sion. In the application to a uniform gravitational field, 
(51) is equivalent to ( [M| ) and (52) is equivalent to (B5). 

The SIVP approach of Sec. |BI D| describes the same 
crossover from rods of size a to rods of size tr -4 1 under a 
repulsive contact interaction of increasing strength. The 
pressure prohle (44) evaluated for h = 0 and any a repro¬ 
duces 


When the same expression is evaluated for 
V = —00 it connects again with (29) but now for cr -4 1 


provided the scaled variables are properly adjusted. 
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IV. POWER-LAW TRAP 

Optical or magnetic traps of several different designs 
for atomic or molecular gases produce wells with a range 
of profiles. How does the pressure at the center of the 
trap vary with temperature? How does the shape of the 
trap potential affect the profiles of density and pressure? 
Power-law traps are well-suited for our two approaches 
and can illuminate these questions with answers from an 
exact analysis. 


A. Profiles for lattice and continuum 


In this application we consider an infinite row of cells 
numbered i = 0, ±1,±2,... at positions Xi = iVc- The 
rods are confined to a region centered at x = 0 by the 
symmetric power-law potential 


u{x) = Uq 



Of > 0 


(53) 


with uo representing a depth and Xq representing (at least 
for a > 1) a width of the trap. 

The analysis proceeds as in Sec. HI A Expression (28) 
remains unchanged. However, the pressure profile is now 
determined by the relation 


exp 


(a — l)(p — pT)\ — I 


crT 


^VtIctT _ Y 




(54) 


where, in addition to p and T from (26), we use the scaled 
variables 


Pt = 


Pt 
P s ' 


(55) 


where fniz) is the Fermi-Dirac function. 

The rods remain confined at any finite temperature: 
we have pr > 0 if T < oo for any a > 0. The limit 
a —>■ oo, representing a trap of width 2xo with rigid walls, 
is subtle. Relation (58) reduces to 


r(I) ^ = 


— = I. 

Xs 


(60) 


We then have p = px inside the trap, with a uniform 
density, pi™®®) = Xg/xg. We thus recover the familiar 
EOS pVc/kx,T = — ln(l — p(™®®)) of the ideal lattice gas. 

The limit T —)■ 0 yields pr = 1 and the profiles for 
pressure and density are 

p= (1-|xr)6»(I-x), pi”®") = 6»(I - x), (61) 


respectively. If a > 1 (a < 1) the pressure px at the 
center of the trap increases (decreases) with T rising from 
zero. For the linear potential (a = 1), which is equivalent 
to the case of the uniform gravitational field (Sec. HI A), 
we have px = ^ for all T. 

In Fig. [^we show pressure and density profiles for the 
two extreme rod sizes in two different power-law poten¬ 
tials. The shape of one potential (a = ^) is concave and 
that of the other {a = 2) convex. Corresponding profiles 
for a linear potential (a = 1) have already been shown 
in Fig. [^albeit on a somewhat different scale. 

The opposite T-dependence of the pressure near the 
center of the trap is evident. The curves at th e lowest 
temperature are close to the T = 0 profiles (61). Natu¬ 


rally, the pressure decrease with rising T in for a = | has 
a much larger effect on the density than does the pres¬ 
sure increase for a = 2. Another noteworthy feature is 
that the difference in profile between the lattice system 
and the continuum system is much more pronounced for 
the convex potential than for the concave potential. 


A solid stack of rods extends out to Xg = ^NaVc and the 
pressure at the center becomes Ps = uo{xs/xo)°‘/{aVc). 
In the continuum limit, a ^ oo, 14 —>■ 0 with cr 14 = 14, 
Eqs. (28) and (54) turn into 


P 


(mes) _ 


1 


l + T/p 


Pt — P — Tin ^ = |x|“. (56) 

P 


One additional relation is needed to bring closure to 


(28) and (54), namely 



dxp)”®")(x) = 1, 


(57) 


reflecting mass conservation. The pressure px x = 0 
can be determined from this relation. For the two ex¬ 
treme rod sizes we thus obtain 


r(l/a + l)/i/„(eP-/^ - l)f = 1 (a = 1), (58) 

= 1 (cr = oo), (59) 


Px 



f 

/ dk 

1 — k — — In A: 

Jo 

Pt 



FIG. 7: (Color online) Profiles of pressure (left) and density 
(right) in a power-law trap with a = 0.5 (top) and a = 2 
(bottom) at temperatures T = 0.02, 0.5,1 for rods of sizes 
CT = 1 (solid curves) and cr = oo (dashed curves). 

































10 


B. Oscillations in density profile 


In Sec. Ill B we have already identified some hard-floor 
effects in the form of spatially attenuated oscillations in 
the density profiles and pi of rods and mass, respec¬ 
tively. We found (in Fig. that the effect is very con¬ 
spicuous in the former and partially averaged out in the 
latter. A more systematic analysis of hard-wall effects 
for rods of various sizes on the lattice and for rods in a 
continuum will be presented in Sec. |Vj 

Here we briefly examine the question whether the soft 
walls of power-law traps also produce patterns of spatial 
density oscillations. We begin with rods of size tr = 2. 
In Fig. we compare data for traps with walls of two 
different degrees of softness. These data are produced by 
using the pot ential ( [5^ with x = iVc and calculating the 
profiles from (BIO) and (10). We use M = 2Jniax + l with 
Imax = 20 and a shift that positions the rod with index 
I = 0 at the center of the trap. 

The harmonic trap is, in some sense, the smoothest 
form of confinement. We see in panels (a) and (b) that 
oscillations do make their appearance in the hi-proflles. 
The amplitudes of these oscillations tend to be rather 
uniform across the region where the rod population is 
significant. The general trend is that the amplitudes in¬ 
crease with increasing chemical potential p, = p/uo or 
decreasing temperature T = ksT/uQ. 

Superimposed on this systematic trend is an oscillatory 
dependence on p. The data in the panels (a) and (b) are 
for two successive values near maximum and minimum 
amplitude of hi-oscillations. The minimum amplitude is 
almost imperceptible. We also note that in the pi-profiles 
these oscillations are almost completely averaged out. 
Increasing the stiffness of the trap walls produces no 




FIG. 9: (Color online) Density profile fii for rods of size cr = 3 
in a power-law trap with a = 10. We use xo = 2014. All 
data are for T = 0.1. The scaled chemical potentials are (a) 
p = 1.1,1.15,1.2 and (b) p = 2.1,2.14,2.2. In each plot the 
second and third data sets are vertically displaced by 0.2 and 
0.4, respectively. 


dramatic changes. The dependence on p of the spatial os¬ 
cillations remains qualitatively similar. The amplitudes 
still oscillate as p is varied. The data in the panels (c) 
and (d) are for a = 10 and for values of p near successive 
maximum and minimum amplitude. 

One systematic trend as the trap wall becomes increas¬ 
ingly firm is that the spatial oscillations become weaker 
near the center of the trap and stronger near the walls. 
This trend is visible between the data in panels (a) and 
(c), for example. We shall see in Sec. that the oscilla¬ 
tory dependence of the amplitudes on p disappears when 
the walls become hard (a —f oo). 

In Fig. [^we show some data for rods of size cr = 3 in 
a trap with relatively stiff walls (a = 10) and twice the 
width of the o ne us ed be fore. The calculations are based 
on Eqs. (B15) and (B16). Here we only show data of the 
hi-profiles for two conditions. 

The average numbers of rods in the data of panel (a) 
are such that the rods easily fit into the space at low 
potential energy. We observe significant oscillations with 
a period near but not exactly three. A small change in 
chemical potential produces shifts in the oscillatory pat¬ 
terns and variations in the average amplitude. We also 
observe slightly larger amplitudes near the wall compared 
to the center of the trap. At significantly smaller values 
of {N) the oscillations are much weaker. 

The data in panel (b) are for circumstances where the 
rods are squeezed into the trap. Here the rod positions 
are more correlated. The oscillations are closer to period 
three. We also observe the oscillatory dependence of the 
amplitude on p. For some values there are one or two 
dominant configurations, producing high amplitude. For 
other values, there are three configuration that have very 
similar statistical weight, producing low amplitude. 


FIG. 8: (Color online) Density profile of rods fii and mass 
Pi for rods of size ct = 2 in a harmonic trap (a = 2) and in 
a power-law trap with firmer walls {a = 10). In both traps 
we use Xo = 1014. All data are for T — 0.1. The chemical 
potential p has the values (a) 1.0, (b) 1.55, (c) 0.3, (d) 0.5. 


V. STERIC WALL EFFECTS 

Microscopic density profiles of rods or other particles 
with shapes near hard walls are relevant in the contexts of 
granular matter, porous solids, and zeolites among oth- 
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ers. Robledo and Rowlinson [3D] studied the effects of 
confinement on hard rods in a continuum. Davis |9] ex¬ 
tended that study to include first-neighbor interactions 
between rods. More recently, Ibsen et al |3Dj reported a 
general and exact solution for hard rods confined by a 
gravitational field and a hard floor. A computer simula¬ 
tion study of Mehrotra et al EH of hard spheres under 
the same confinement in 3D produces similar results. 

Here we pick up threads from Secs. IIIB and IV B to 
investigate steric wall effects of hard rods on a lattice. 
One goal is to showcase the versatility of the method of 
analysis presented in Appendix [^ and to establish how 
it connects to the continuum analysis familiar from pre¬ 
vious work. 

We consider a box with rigid walls and investigate the 
oscillations in the density profiles produced by the steric 
interactions between rods of size cr > 2. The effects of 
a single wall, relevant in sufficiently wide boxes, can be 
determined analytically for rods of any size on the lat¬ 
tice and for rods in the continuum. Two walls within 
the distance of a certain coherence length affect the den¬ 
sity profile from opposite sides. That coherence length 
is shown to grow with average density. We combine ex¬ 
act analytic results with results from a rigorous recursive 
scheme. 


A. (7 = 2 


We begin with the case of a semi-infinite box with one 
wall at * = J and the other at * —>■ —oo populated by rods 
of size (T = 2 to an average mass density 0 < < 1. 

Later we move the second wall to i = —I. For the semi¬ 
infinite box we have found an analytic solution. The 
result (with j = I — i in the present context) turns out 
to have a simple structure: 



(mes) 


1 - 


p(mes) 
p(mes) _ 2 


J 


j = 0,l,2 ,... 


(62) 

The relation between the fugacity and the average mass 
density is 


p(mes) _ ^ _ 


ymc 


or C =4 


(1 _ p(mes)^2 


- 1 


(63) 

We have derived ni in (62) directly from ( Bllj ) with g 2 in 
the form of an infinite continued fraction that is readily 
evaluated. The hi for i = 2, 3,... then follow directly 
from (B5|. The oscillations in the probability distribu¬ 


tion Tii of rod positions thus decay exponentially with 
distance from the wall. The boundary coherence length, 




1 


ln(2/p(“®'^) - 1)’ 


(64) 


vanishes for 0 and diverges for 1. 

Next we examine how the oscillations near the wall at 
i = I are affected by the presence of a second wall at i = 



FIG. 10: (Color online) Density profile for rods of size (7 = 2 
with center at coordinate i = 0,±1 ,..., / in a box of width 
/ = 10. Full circles con nected by lines: solution of (BlOl. 
Open circles: result (621 for a box of infinite width with one 
wall at i = 10. 


—I. For that purpose we have solved (BIO) with C, from 
(63) for comparison with the analytic solution (62). The 
results, shown in Fig. |10[ are almost indistin guis hable for 
^(mes) < Q_ 75 _ Here the coherence length (64) is much 
smaller than the distance between the walls. At larger 
mass density the two sets of results begin to deviate from 
each other as ^ grows and reaches about half the wall-to- 
wall distance at < 0.95. 


The signature wall effect for ct = 2 manifests itself in 
the form of attenuated, period-2 spatial oscillations in 
the hi. The exponential attenuation is governed by a 
coherence length that grows with the density of rods in 
the box. In the limit p(“®®) —>• 1 the oscillations persist 
across the box as expected. 


B. 2 < (7 < oo 


We continue with the analysis of a system of rods of 
size 2 < (T < oo subject to the potential 


U, = 


oo : i < 0, 

0 : z > 0, 


(65) 


representing a single rigid wall at z = 0. We present the 
exact solution of Eqs. (B2) adapted to this case. 


2+cr—2 

k 

i+fT—1 

k 

C = hi 

1 - II ^3 

n 

1 - H ^3 

k—i 

j=k—iy+2 

k—i 

j=k-cr+l 


(66) 
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We have = 0 for i < 0 and assume that far from the 
wall (i ^ 1) the hi approach uniformity, 


lim hi = ^ 


(mes) 


(67) 


where 0 < < 1 is the mass density in the bulk 

(average cell occupancy). The fugacity, the only control 
variable aside from a, depends on the asymptotic solution 
(1^ as follows: 


c = 


fj(mes) [1 (cr - l)h(“®®)] 


<7 — 1 


[1- 


an 


S)]' 


0 < 


( 68 ) 

Next we convert ( |66[ ) into a recursion relation that ex¬ 
presses the solution at a given site as a function of the 
solutions at the cr — 1 sites immediately to its left: 

i-l 

hi = 1 - hj 

j—i—G-\-l 


CT+l 

C 


i-l 

n 

k—i—G+l 


1- E 


j—k — G-l-2 


1- E ' 

j=k — /7+l 


(69) 


for i = a, a + 1,... and with ( from ( 681 . This recursion 


relation depends on the cr — 1 parameters ni,..., ncr-i- 
Assuming that the solution with asymptotics (67) is 


unique, these parameters can be found with a little guid¬ 
ance from (|62| for the case cr = 2. They are 


hi = - 


[ 1 - 


^(mes)jr 1 


an 


[1 - (cr - l)h(“®®)]’ 


i = (70) 


In Fig. [D we show data generated recursively from 


(69) for the scaled density i/i = arii versus the scaled 


position Xi = {i — l)/cr for rods of sizes cr = 5 and 
cr = 200. The latter size is meant to generate an im¬ 
pression of what to expect in the continuum limit. The 
data suggest that the dominant wall effect manifests it¬ 
self again in an attenuated spatial oscillation. Succes¬ 
sive minima are approximately spaced by cr. Only in the 


limit —)• 1, when the attenuation disappears, are 

the oscillations locked into the wavelength cr. This raises 
the interesting question, best analyzed in the continuum 
limit, what the spectrum of the wall oscillations is and 
how it depends on 


C. cr = (X) 


The continuum limit carried out for 0 < a; < 1 pro¬ 
duces an exponential function as follows: 


v(x) = se s = —-—-r. 

^ ^ ’ 1 - p(mes) 


^(nies) 


(71) 


This result can now be extended to x > 1 by using a 
continuum version of the fugacity (68), 


C = se®. 


(72) 


and a continuum version of (66) 


h(x) = Cexp^—y dx'h{x')^ , (73) 


where 


h{x) = 


iy{x) 


1 - dx' v{x') ‘ 


(74) 


as derived Percus Following Vanderlick et al. [33], 
we convert (73) and (74) into difference-differential equa¬ 
tion. 


dx 


h{x) = h{x) \h{x — 1) — h{x )\, 


(75) 


dx 


v{x) 


h{x) 


= iy(x) — v{x + 1), 


(76) 


respectively. Next we integrate (75) and (l7|) in alter¬ 
nating sequence over intervals of unit length, using (71) 


and h[x) = 0 for x < 0. We thus obtain the following 
exact, continuous, and piecewise analytic expression for 
P on successive intervals m < x < m 1: 



FIG. 11: (Color online) Scaled density of rods of size ct = 5 
(big circles) and a = 200 (small circles) with first cell at 
position Xi = {i — l)/o' near a rigid wall at i = 0. 


"L k+i 

= (77) 

k^O 

with asymptotic value, 

lim il(x) =-^ (78) 

x^oo 1 + S 

far from the wall, approached more and more slowly with 
increasing average mass density The lattice DFT 

analysis thus connects neatly with known results [9||29l 
|30l |32l |33| for monodisperse rods in a continuum. 

In Fig.[^we graphically compare this analytic solution 
([T^ for the continuum model with the iterative solution 
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FIG. 12: (Color online) Scaled density of rods of size ct = 20 
on a lattice (big circles) and rods of (scaled) unit size in a 
continuum (solid lines). 


(69) for the lattice model with a = 20. The dominant 


feature of the curves is an attenuated spatial oscillation. 
The singularities at a; = m become progressively weaker 
with increasing m: d^v{x)/dx^ is continuous at a; = m 
for k < m. Except for to = 1, the singularities do not 
coincide with the minima of i'(x). 

The limit —)■ 1 (s —>■ oo) is subtle. As the popu¬ 

lation of rods becomes more crowded, they begin to line 
up with increasing p rob ability near the integer positions 
as illustrated in Fig.[l^a). In the limit 1, max¬ 

ima and minima approach the singularity values x = m 
in pairs from opposite sides. The maxima diverge and 
the minima approach zero. The area under the curve 
between successive singularities approaches unity. The 
shape of the curve approaches an L of infinite height and 
unit width, effectively the function 


lim ^{x) = y 

s—>-oo 

m—Q 


5(x — to). 


(79) 


In Fig, p^b ) we show the spectrum of the attenuated 
oscillation (|77|), 


P{q) = 2 


dx 


D{x) - COs{qx). 


(80) 


This quantity features a peak at wave number 
1 < g/TT < 2 emerging from a broad and flat background. 
The Fourier transform (801 and the limit s —>■ oo are 



FIG. 13: (a) Scaled densit y o f rods in a continuum with 
0.9. (b) Spectrum (80) of the attenuated oscillation 


(mes) _ 


(771 at three values of 


not interchangeable operations. With increasing p(“®®) 
the peak becomes taller and sharper as it moves toward 
commensurability at qj-K = 2. 

The DFT mass density in the continuum, calculated 
from via a continuum version of (|10[), 


p{x) = / dx'v{x')^ 


(81) 


becomes 


p(^) = 1 - 


k—Q 


k ^—{x — k)s 


(82) 


on successive intervals TO<a;<TO-|-l. In Fig. |14| we 
show profiles of p{x) for the corresponding to the profiles 
of v{x) shown in Figs. 12 and 13 The oscillatory wall 
effects are somewhat milder in p{x) than in i>(x) but still 
fairly conspicuous. In the limit s —>■ oo the oscillations 
in v{x) diverge as in (79) but disappear in p{x), which 
becomes constant. 



FIG. 14: Profiles of mass density of rods in a continuum at 
three different average values 


VI. CONCLUSIONS AND OUTLOOK 

In this work we have investigated the equilibrium sta¬ 
tistical mechanics of monodisperse hard rods confined by 
external potentials, producing characteristic spatial pro¬ 
files of density and pressure. The rods occupy a consec¬ 
utive cells of volume 14 in a linear array. The continuum 
limit is implemented as u —f oo, 14 0 with aVc kept 

finite and nonzero. The rods interact via hardcore exclu¬ 
sion forces. Contact forces have also been considered. 

We have been employing two distinct methods of anal¬ 
ysis, density functionals and generalized exclusion statis¬ 
tics, with partially overlapping domains of applicability 
and somewhat complementary strengths. The usefulness, 
soundness, and consistency of the two approaches has 
been demonstrated in a series of applications that in¬ 
clude confinement by rigid walls, gravity, and power-law 
traps. 

We have shown that profiles of density, for example, ex¬ 
hibit characteristic features on a mesoscopic length scale 
that are reproduced identically by both approaches and, 
for tr > 2, additional features on a microscopic length 
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scale that are only resolved in the DFT analysis. The 
DFT pressure profiles are also compared with profiles of 
average microscopic pressure inferred from pair distribu¬ 
tion functions in the context of one application. Finally, 
we have established contact between our results for the 
lattice system and prior results for corresponding contin¬ 
uum systems, demonstrating full consistency 


becomes 


PmUk, t)=Y, (A2) 

TJlr 

n 

' / ^ J mn\^/ ■, /,\ ^ 

m KXmn{t) 


The statistical mechanical analysis of hard rods of 
mixed sizes presents itself as a natural and promising ex¬ 
tension of this work. Significant ground has already been 
broken via both calculational tools used in this work. 
The mathematical structure of exact density functionals 
for hard-rod mixtures has recently been established [15] 
and is awaiting analysis in specific applications. 

Hard rod ensembles of mixed sizes on a lattice are 
mathematically equivalent to magnetic domains in Ising 
chains. The combinatorics and statistical mechanics of 
such domains analyzed via generalized exclusion statis¬ 
tics have been in place for some time |34j and is awaiting 
extensions to heterogeneous environments and new appli¬ 
cations of current interest. The completion of the work 
reported here has paved the way for projects along these 
lines [35]. 


To obtain the corresponding expression in real space, we 
calculate 


dfc 


— ikXmri-\-'ikx 


sin(fca 


kxr, 


dk 


— ikX^ji-\-ikx 


1 ikXmrjO^ 

da e 


da(5(x X-uiYi CXX-uin) 


2|a 


' 5(^7 ^n) j 


(A3) 


where 


g{x;xm,x„) 


1, a; e [min(a;„,a;„),max(a;m,a;„)], 

0, else. 

(A4) 


Appendix A: Local pressure from momentum flow 


In terms of the particle (hard rod) momenta n„ (t) and 
positions Xnif), the microscopic momentum density in 
Eq. ^ is given by n(a;,t) = J2n^ri{t)S{x - x„(t)) . 
The time derivative of its Fourier transform n(fc,t) = 
/ da; n(a;, t) nn(i) reads (with rur 

the mass of the rods, and assuming only pair interaction 
forces to be present) 


Accordingly, 


Pmic {x,t) = 'y' S(^x-Xn (t) ) (A5) 

TTlr 

n 

+ I frnn{t) I 9{X; Xm{t),Xn{t)) . 


Using fmn ^mn — fnm^nm 
can be written as ^ 


and Eq. (]A4[), the double sum 
and the 


x^<x ^x„>x fmn-xmn 


Xr 


equilibrium average of Eq. (A5) then yields Eq. (231. 


5n 



g-ikx^{t) 


(Al) 


n,m 


+E/“‘w 


^-ikx„(t) 


Here fmnit) = /(a^m(f), a^n(i)) is the force of the mth on 
the nth particle (/„„ = 0), /®’"*(t) = /®’"*(x„(t)) 
the external force on the nth particle, and we 
have used in = H^/mr and Newton’s equations 
Hn = /r‘ +E^/n™- For interaction forces obeying 
the principle of action and reaction, fmn = —f nrm the 
second term on the right hand side of Eq. (Al) can be 
written in the form (1/2) Y.n,m fmn (e"*''"" - = 

-ikY,m,nfmnXnine.~'‘^^^'^[s\n{kXnin)/{kXmn)\, where 
^mn = ( xrn+Xn)/‘^ and Xmn = (x„ —Xm)/2. Comparing 
Eq. ( |Al[) w ith the Fourier-transformed right hand side 
of Eq. ([22[), the microscopic pressure in Fourier space 


Appendix B: Density profiles 

Consider a system of rods of size a confined to a fi¬ 
nite array of cells numbered l,...,L = M-|-cr—1. The 
hardcore repulsion is the only interaction between rods. 
The external potential Hi is arbitrary. The free energy 
functional @ simplifies into 

L ( i i 

/3E =^“1 -I-^1- E E 

i—1 \ j—i—G+1 j—i—a+l 

2-1 2-1 'I 

j=i — U + l j=i — (y+l ) 

where fii,, um are the (yet undetermined) probabili¬ 
ties of rods at the allowed positions. The minimization 
© then leads to the following relations that determine 
the exact density profile of rods, {fii}, for an any given 
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external potential lAi at given temperature T and chem¬ 
ical potential fj,: 




2 +cr—2 

k 

n 

1 - ^ fij 

k—i 

j=k—<7+2 

■i+cr—1 

k 

n 

1 - ^ fij 

k—i 

j=k—<7+1 


(B2) 


for i = 1,... ,M with hi = 0 for i < 1 and i > M im¬ 
plied. In some applications we use a semi-infinite array 
(M —> oo) with the second boundary condition replaced 
by a prescribed (zero or nonzero) limit hoc- In the fol¬ 
lowing we use the control variables, 




(B3) 


with the implication that Ai = 0 for i < 1 and i > M. 

For rods of size cr = 1 the hardcore repulsion does not 
produce any interference between the hi at different po¬ 
sitions. Equations (B2) remain uncoupled. The density 


profile rii of rods, which, in this case, coincides with the 
mass density profile pi, reads 

h, = -4^: z = l,...,M. (B4) 

Its only structure is that imposed by the potential Ui. 

In the case cr = 2 a rod at position i obstructs the 
placement of a rod positions i — 1 and i + \. This inter¬ 


ference is reflected in Eqs. (|B2|), which now read 


i{l - ni) 


{l-m- -rii- Ui+i) ’ 


f = 1,..., M, 

(B5) 

with np = riM+i = 0 implied. We solve the coupled 
Eqs. (B5) by a strategy that also works for cr > 2 as we 


shall see. We reduce the set of nonlinear equations into 
two sets of recursion relations to be solved in sequence. 
This method has the benefit of isolating the physical so¬ 
lution. We begin by introducing the auxiliary quantities. 


hi = 


1 Tii 


i = l,...,M. 


(B6) 


We thus convert (B5) into 


= K + hihi-i, i = l,...,M, (B7) 

from which we determine the hi recursively: 

CA, 


hi = 


1 -|- hi- 




i-l 


(B8) 


The rii are then, in turn, determined recursively via (B61: 

h 


rii = 


1 + hi 

producing the explicit form 


(1 - n^+i), z = M,..., 1, (B9) 


M-i 


I 


1)^ n 1 


h 


i+k 


-, i = l,...,M. (BIO) 


1=0 k=0 

A different rendition of that solution is derived from 


(B9) with use of (B8): 


rii = 


hi 


i = 1,...,M, 


1 hi gi+i 

where the gi are generated recursively: 

CK 


(Bll) 


9^ = 


1 + gi+i 


, i = M,...,l. 


(B12) 


If the external potential is symmetric under under re¬ 
flection, as in power-law traps or boxes with rigid walls, 
we have Xm+i-i = Ai, i = 1,... ,M. The density pro¬ 
file of rods must then also exhibit that symmetry. To 
make this symmetry transparent we recognize that we 
have gM+i-i = hi, i = under these circum¬ 

stances. We can then transform fiM+i-i into hi as fol¬ 
lows: 


iM+l-i 


(Xm+I-i 

1 + hjif-i 


CAm+i-i 


1 + hM+l-i + gM+2-i 1 I C^M+l-i {1 + gM+2-i)i^ + hM-i) + C^M+l-i 

t + . . , -h gM+2-i 


1 -|- hM-i 

CA* 


CA* 


1 -|- hi i 


(1-I-/ii_i)(l-I-5i+i)-I-CAi 1 I „ I ‘’Ai 

-L + gi+i H- 


1 + ^i + gi+l 


1 + hi-i 


(B13) 


Now we generalize this method to rods of unrestricted size on the lattice: cr = 1,2,.... A rod of size cr at 
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position i obstructs the placement of rods at positions 


i±l,...,i±(cr—1). The right-hand side of (B2) contains 
hi at 2(7—1 consecutive positions. The coupled equations 
now involve polynomials of order a. The boundary condi¬ 
tions require that we set hi = 0 for i = —cr -1- 2,0 and 
i = M-|-l,...,M-|-CT—lin these M coupled equations. 
The auxiliary quantities 


hi = 


<7— 1 


1 ^ ^ ‘^i+k 


-1 




then convert (B2l into the set of recursion relations 

(7 — 1 

hi=CX^Y\{^ + hi-k) \ (B15) 




For given hi the rii then follow recursively from (B14): 

T / (7 — 1 \ 


1 + hi 


1 ^ ^ 

V fe=l / 


i = M,...,l. (B16) 


Appendix C: Pair distribution function 

Consider a system of N hard rods of size a at positions 
i/c, k = 1,... ,N on a ID lattice of L sites. The system 
is confined by hard walls at positions io and iN+i- The 
lattice partition function 

2iV + l—<7 2iV + l—2tT 

^Jv(*oDAr-l-l) = ^ ^ 

iN—'i'0+N(7 iN_i—iQ-\-{N — l)(7 
iN+i-Na 

X (Cl) 

ii=io+<y 

with variable change jk = ik — io — ka and excess volume 
defined as 

Lex = L - Na = iN+i - io - {N + l)a, (C2) 


becomes 

Ziv(zoDiv+i) = g Y 


Lpx Lp 


X Y^ A 

ii=o 


liv=OiN-i=o 

(C3) 


(^14) and, after summation over jk, 


I _ g-/3mg(Lei.-|-l) 
1 — e-Prng 


N 


(C4) 


The density and pair distribution functions are then ex¬ 
pressible as follows in terms of partial partition functions 
Zk{l,rn) representing systems of k rods confined by hard 
walls at positions I and m min Haul]: 


p{i) = 




^-/3mgi 


iN+l) 


^ ^fc-l(*OD)^Af-fc(b *Ar-l-l)j (C5) 


k=l 


g-pmg{i+j) 


iN+l) 

N 


(C6) 


X Yi ^k-i{io:i)Zi-k{i,j)ZN-i{j,iN+i)- 


k<.l — l 


A lattice version of expression (24) for the average mi¬ 


croscopic pressure to be used in Sec. Ill C for comparison 


with the DFT pressure profiles thus reads 

(7 

P-micii) = ksTp{i) + ksTY^ P^‘^\i - h,i + o - k). (C7) 


k=0 
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The averages are over the (grand-canonical) Boltzmann 
distribution for Hamiltonian ||^, but with Lli replaced by 
^Mermin^ whcrB is thc unlque external ‘Mermin 

potential’ that (for given interaction Vij) would generate 
the hi as the mean occupation numbers in equilibrium. 
As shown in Ref. m the Oi, bij and dij refer to the prob¬ 
abilities to have zero or one rod in ‘one-particle cavities’, 
i.e. in intervals, where at most one rod can be placed. 

In one dimension, pressure has the same dimension as 
forc e, namely energy per length. The scaling factor 14 in 
(131 makes p{x) dimensionless, consistent with (12l. 

In the case presence of phase transitions only the latter 
combination may be unique. 

The microscopic pressure at £ = 0 is not equal the total 
weight of the rods but it has the form po = -(-pn, 

as can be shown from (C4l, where pl is the pressure at 
the top of the system. With increasing system size pl 
approaches zero. 





